DESY 10-061 

SHEP-10-04 

CERN-PH-TH/2010-091 



Using HERA Data to Determine the Infrared 

Behaviour of the BFKL Amplitude 

(N 
> 

j^ : H. Kowalski \ L.N. Lipatov ^'3, D.A. Ross \ and G. Watt ^ 

^ I ^ Deutsches Elektronen- Synchrotron DESY, D-22607 Hamburg, Germany 

^ Petersburg Nuclear Physics Institute, Gatchina 188300, St. Petersburg, Russia 

^ II. Institut fiir Theoretische Physik, Universitdt Hamburg, D-22761 Hamburg, Germany 

^ School of Physics and Astronomy, University of Southampton, 



X 



5 



Theory Group, Physics Department, CERN, CH-1211 Geneva 23, Switzerland 



^\ 

«D : Highfield, Southampton SO 17 IB J, UK 

^, 

(N 
> 

in 
vn 

m 

p 

o 
o 



Abstract 

We determine the infrared behaviour of the BFKL forward amphtude for gluon-gluon scatter- 
ing. Our approach, based on the discrete pomeron solution, leads to an excellent description 
of the new combined inclusive HERA data at low values of x (< 0.01) and at the same 
time determines the unintegrated gluon density inside the proton, for squared transverse 
momenta of the gluon less than 100 GeV^. The phases of this amplitude are sensitive to the 
non-perturbative gluonic dynamics and could be sensitive to the presence of Beyond-the- 
Standard-Model particles at very high energies. 
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1 Introduction 

One of the major results from HERA (see [1], and references therein) is that the inclusive 
cross-section for the scattering of virtual photons against protons at low x (i.e. high energy), is 
dominated by the gluon density inside the proton. This allows one to study the behaviour of 
the gluon density as a function of gluon momenta, i.e. the fraction of the proton's longitudinal 
momentum x and the transverse momentum k. The study of the dynamics of the gluon 
density is usually motivated by its importance to other physics reactions, like dijet or Higgs 
production at the LHC In addition to this merely "utilitarian" aspect the dynamics are very 
interesting because the gluon density is a fundamental quantity, comparable to black-body 
radiation in QED, and because gluon-gluon interactions are the source of the forces which 
keep matter together. 

The dynamics of the gluon distribution at sufficiently low x is best determined by the 
amplitude for the scattering of a gluon on a gluon, described by the BFKL analysis. In this 
analysis the pomeron is considered as a composite state of two so-called reggeized gluons [2] . 
One of the salient features of the purely-perturbative BFKL analysis is the prediction of a 
cut-singularity with a branch point A leading to a low- a; behaviour for the gluon density of 
the form 

xg{x) ~ x~'^, (1.1) 

with only logarithmic corrections in x. In leading order, A is given by 

121n2 , , 

A = as. (1.2) 

TT 

The branch point A only depends on Q^ through the running coupling a^. Experimentally, 
this branch point is given by the rate of rise of F2 with diminishing x, F2 ~ x^^. Thus, for 
many years, it was claimed that BFKL analysis was not applicable to HERA data, firstly 
because the value of A obtained from (II. 2p was much larger than the observed value, and 
secondly because HERA found substantial variation of A with Q^. 

The first of these difficulties was ameliorated by the NLO contribution to A |3] , once the 
very large corrections were resummed using the collinear resummation technique [1]. On 
the other hand, the second difficulty, namely the question of the large Q^ dependence of the 
parameter A remains problematic. In this paper, we will show that using a modification of 
the BFKL formalism, which leads to discrete solutions (i.e. Regge poles rather than a cut), 
the Q^ dependence of A can indeed be reproduced from the BFKL amplitude. At the same 
time an excellent description of the low-x HERA data is obtained, thereby determining the 
infrared behaviour of the BFKL amplitude. 



2 BFKL analysis 

The fundamental ingredient of the BFKL analysis is the amplitude for the scattering of a 
gluon with transverse momentum k off another gluon with transverse momentum k' at centre- 
of-mass energy i/s which is much larger than the momentum transfer and much larger than 



the magnitudes of the gluon transverse momenta. In the forward case (zero momentum 
transfer) this amphtude, A{s, k, k'), is found to obey an evolution equation in s given by 
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Ais, k, k') = 6ik^ - k'^) 6 (in ^)+ f cig'/C(k, q)^(s, q, k'), (2.1) 



where /C(k, k') is the BFKL kernel, currently calculated to order a^ [3]. The kernel is obtained 
by summing all graphs which contribute to this process, but keeping only leading (and sub- 
leading) terms in Ins. Such graphs can be drawn in terms of an effective "gluon ladder". 
The Green function evolution equation (12. ip can be solved in terms of the eigenf unctions of 
the kernel 

j dk''}C{k,k')fUk') = cufM- (2.2) 

In leading order and with fixed strong coupling a^ the eigenf unctions are parameterized by 
a "frequency" u and are of the form 

/.(k) = {kT''^\ (2.3) 

with an eigenvalue, oj, given by 

i^ = asXo{i^), (2.4) 

where Xo{^) is the leading order characteristic function. The maximum value of u, ai u = 0, 
is equal to the branch point, A. 

In a recent paper [S], we reported on an attempt to fit HERA data for structure functions 
at low values of Bjorken-x, using the first few discrete solutions of a modified BFKL [2J 
equation. The modification was proposed by one of us [6j^ in order to obtain discrete poles 
(as opposed to the above-mentioned cut) thereby attempting to explain how QCD could, at 
least in principle, reproduce the very successful results of Regge theory in hadronic processes, 
and was the progenitor of modern string theory. 

In the modified BFKL approach the strong coupling constant a^ is running as one moves 
away from the top or bottom of the gluon ladder. This allows the transverse momenta of the 
gluons k, which dominate the amplitude, to have a large range as one moves away from the 
ends of the laddecl and results in a solution to the eigenvalue equation, (12. 2p . in which the 
frequencies of oscillation u are themselves k dependent. In the semi-classical approximation, 
valid in the region 

^ In(Kk)) 



dink 
the solution for a given eigenvalue u has the form 



« z/(k), (2.5) 



/^(k) = ^li^^expf22 /° u{k')dlnk' ] . (2.6) 



^Discrete solutions in order to obtain discrete poles of the BFKL equation were first discussed in [21[71[S]. 
^Here the transverse momenta are controlled by the convolutions with impact factors that determine the 
precise process under consideration. 



Here the function i^(k) is determined from the perturbative expansion of the kernel /C, 

a,(A;2)Xo(^(k)) + a,2(A;2)XiMk)) = a;, (2.7) 

where Xo ^^^ Xi ^''^^ the LO and NLO characteristic functions, respectively. We note here 
that in our numerical analysis, we modify Xi following the method of Salam [3], in which the 
collinear contributions are resummed, leaving a remnant which is accessible to a perturbative 
analysis. Since the RHS of eq. (12. 7p is constant it follows that the frequency z/ which is the 
argument of these characteristic functions has to depend on k in order to compensate for the 
k dependence of the running coupling. Eq. (12. 7p provides a match to the DGLAP equation [9] 
which, in the limit of small u, has an anomalous dimension 7aj(k) related to i/(k) bjo 

7^ = ^+Mk) + ^. (2.8) 

For sufficiently large k, and for positive values of co, eq. (12.71) no longer has a real solution 
for z/ and the eigenf unctions become exponentially decreasing as opposed to oscillatory. The 
transition from the real to imaginary values of z/(k) singles out a special value oi k = kcriti^), 
such that h'uikcrit) = 0. The solutions below and above this critical momentum fccrit, niust 
be carefully matched. Indeed, at /c = /ccrit, eq. (12.61) is invalid because the coefficient C(z/(k)) 
diverges and an Airy function is used to interpolate smoothly between the two regions k <^ 
/Ccrit and k ^ /Ccrit ^3\- The argument of the Airy function is a function of the generalized 
phase, 

^W^ik')\. (2.9) 

The eigenfunction is given by 

kUk) = Uk) = Ai (^-(^0,(fc))f^ . (2.10) 

This matching of the solution to the BFKL equation in the regions k < k^a and k > fccrit 
means that the BFKL equation contains more information than the DGLAP equation near 
u = 0. Indeed, although the BFKL equation is, in principle, an integral equation, whose 
solution must therefore be sampled over the entire range of k', the quasi-local nature of the 
kernel /C(/c, k'), i.e. the fact that the kernel only has non-negligible support where k and k' 
are of the same order of magnitude, means that the kernel can be written in the form 



^(k, k') = ^ E ^"'^^"^ (ln(kVk'^)) , (2.11) 

n=0 

where the coefficients c„ are given by 

c„=/ rffc'2/C(k,k')^^(ln(kVk'2))" (2.12) 

Jo k'nl"- 

•^Note that the last term on the RHS of eq. (|2.8p arises from a conversion from an s-dependence determined 
in the BFKL approach to Bjorken-x, used in the DGLAP approach. 



rapidly falling at large n. The coefficients Cn are, of course, functions of q;s(/c^), with a 
perturbative expansion in q;s(/c^). This leads to the BFKL equation in the pseudo-differential 
form 

dk"lC{k,k')Uk') = ^cj^— — j Uk) = ioUk). (2.13) 

This means that the integral BFKL equation is equivalent to a quasi-local equation which 
can be cast in the form of the pseudo-differential equation 

k j dk'^1C{kM)U\^) = x(~tj^,asie)]uk) = uflik). (2.14) 

In the semi-classical approximation for which fc-dependence on fu){k) is such thalQ 

' '^'uk)^m('-^^]", (2.15) 



d\n{k) J \ d\nk 

eq. (12.141) looks like the non-linear differential equation 

x[-^^^^^.<k)j = u. (2.16) 

This is equivalent to the DGLAP equation written in the usual linear form 

1^ = ^lyUc^s{k'))Uk), (2.17) 

where I'uj {c(s{k'^)) is the solution to 

xK(as(A;')),«s(A;')) = CO- (2.18) 

Near k = fccrit, where z/ is small, only the second derivative term in eq. (I2.13P is important, 
and the differential equation then has the form of a Schrodinger equation. Moreover, since in 
that region the difference between as^k"^) and as(/c^j.jj is approximately linear in ln(A;^//i;^j,jJ, 
the "potential" of this Schrodinger equation is linear and the solution is an Airy function. 
Away from this region the Airy function is oscillatory for k -C /ccrit and exponential for 
k ^ /ccrit- The oscillatory behaviour of the Airy functions is the same as of the linear 
combinations of the solutions of the DGLAP equation (I2.17p . 

/,(fc)=exp(±#,(A;)), (2.19) 

where i(pu){k) is the generalized phase of eq. (12.91) . This means that the solution of eq. (12.101) 
provides a single interpolating function which reproduces, to a very good approximation, 
the eigenf unctions of the kernel for all values of transverse momentum, k, with oscillation 
frequency given by eq. (12.181) . 

The DGLAP equation is a linear equation which has two oscillatory solutions in the region 
k ^ fccrit- The BFKL equation, eq. (12.141) . can be considered as a quantized version of the 



*This condition is, in fact, equivalent to the condition of eq. (|2.5 



DGLAP equation f l2.16p . The matching of the BFKL solution in the region k ~ fccrit (which 
is outside the region where the semi-classical approximation is valid since the logarithm of 
the oscillation frequency, ln(z/), becomes infinite at kcnt) imposes unique coefficients for the 
two solutions, eq. fl2.19p . which cannot be obtained from the DGLAP equation alone. 

Moreover, the quasi-local property of the BFKL kernel, which allows it to be recast into 
a differential equation (albeit of infinite order), permits the continuation of the matched 
solution from large transverse momenta k < kait down to a value k ^ Aqcd at which point 
the perturbative analysis becomes invalid. Nevertheless, the facility to continue the solution 
near the infrared region means that the unknown effects from the non-perturbative sector 
of QCD can be encoded into information about the phase of the oscillations at some small 
transverse momentum. Thus, for sufficiently low transverse momentum, k <^ /ccrit, this Airy 
function solution has the asymptotic behaviour 

k Uk) r^ sin (^(P^k) + ^y (2.20) 

For small values of k the non-perturbative effects become important and these determine 
the phase of the oscillations at the point k = kQ ^ Aqcd, denoted by rj. In Ref. [6] it was 
assumed that this phase is common to all eigenfunctions. This gives rise to two boundaries, 
one at k = kcnt and the other at k = kQ, and leads, together with the perturbative result of 
eq. ( 12. 9p . to the quantization condition 

Mko)=(n-^J7r + r]7r. (2.21) 

Consistency with both boundary conditions can only be achieved for a discrete set of 
eigenvalues u and leads to a discrete set of solutions. Figure [1] shows as an example the 
values of cOn for the eigenfunctions fui„, with n = 1,2,..., 120 determined assuming that 
f] = 0, which were found by solving the equation f l2.18p numerically. The eigenvalues a;„ are 
approximately related to the eigenfunction number n by a simple function 

^ 

Un ^ — . (2.22) 

l + 0.95n ^ ^ 



Figure [2] shows the values of the logarithms of critical momenta for the same eigenfunc- 
tions. Note that the critical momenta are growing with increasing n very fast, because the 
eigenfunctions oscillate, approximately, as a function of log(A;). Figure [3] shows the first eight 
eigenfunctions, f^^^, n = 1 . . . 8 determined at r^ = 0. 

The quasi-local property of the BFKL kernel serves a further purpose. The scattering 
amplitude for two gluons is obtained from the eigenfunctions of the BFKL kernel in the 
multi-Regge regime, i.e. Ins ^ Infcj, where kj is the transverse momentum of any of the 
gluons exchanged in the t-channel. The transverse momenta of the external gluons are 
limited by the impact factors of the proton and photon, so that the contribution to structure 
functions from gluons with transverse momentum much larger than the photon virtuality, 
Q, is highly suppressed. In the central-rapidity region, it is possible that gluons can diffuse 
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Figure 1: The eigenvalues Un for the eigenfunctions f^^, with n = 1,2, 
indicates the approximate relation, Un = 0.5/(1 + 0.95 ri). 



120. A thin line 



into gluons with much larger transverse momenta, but in this case the quasi-local property 
of the kernel leads to a suppression of the contribution from such gluons. Thus we are always 
forced to remain within the multi-Regge kinematic regime. On the other hand, since the 
BFKL amplitude is valid for any energy, at sufficiently large energies, since the values of fccrit 
are very large, there would be non-negligible contributions from gluons with large transverse 
momenta, above the threshold of any possible physics Beyond the Standard Model (BSM). 

In [S], we expressed the low-x structure function of the proton, ^2(0;, Q^), in terms of the 
discrete BFKL eigenfunctions by 



F2(x,Q^) = / dz —^Y)is{z,Q,k)xg(^-,k 



(2.23) 



where xg (f , fc) denotes the unintegrated gluon density 



xg{x,k) = ^ 



f*.(*') 



k' X 



k'f:(k')f^M 



(2.24) 



and $p(fc) denotes the impact factor that describes how the proton couples to the BFKL 
amplitudes at zero momentum transfer. The impact factor, ^Y)is{z,Q,k), which describes 
the coupling of the virtual photon to the eigenfunctions is given in |10]; the dependence on z 
reflects the fact that beyond the leading-logarithm approximation, the longitudinal momen- 
tum fraction, x of the gluon differs from the Bjorken-value, determined by Q^. $dis(z, Q, k) 
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Figure 2: Logarithms of the critical momenta /Ccrit for the eigenf unctions f^]^, with n = 
1,2,. ..,120. 

of Ref. [lO] is determined taking into account kinematical constraints allowing for non-zero 
quark masses. 

We were able to demonstrate that a good fit to data could be obtained using only the first 
four discrete eigenfunctions. In that paper [5], however, no attention was paid to the exact 
form of the proton impact factor. In the fit, the unintegrated gluon density was expressed 
as 

xg{x,k) = 5^a„a;--"A;(2+-")/^„(A;). (2.25) 

n 

and only the four coefficients a„ were fitted to the data. The proton impact factor resulting 
from the coefficients a„ turns out to be highly unstable and becomes negative for sufficiently 
large values of transverse momentum. 

In this paper, therefore, we fix the proton impact factor $p(k) in the form 

%{k) = Aee-"'\ (2.26) 

which vanishes at small k"^ as a consequence of colour transparency, is everywhere positive 
and decreases for large transverse momentum, with a maximum to be determined by the 
fit. The particular functional form of the impact factor is not very important as long as it 
is positive and concentrated at the values of fc < 0{1) GeV. The coefficients a.„ are now 
determined from the impact factor through 



«" = E/f^*^(^')(^')'""/i(^')- 



(2.27) 



In order to do this, it turns out that it is necessary to take many more eigenfunctions 
because the values of the overlap integrals in eq. 02.27P are diminishing very slowly, a„ ~ 
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Figure 3: The first eight eigenfunctions /^ 



n 



, 8 determined at r^ = 0. 



\l \pn. This is due to the fact that the amphtude of the eigenfunctions is dropping hke 
~ 1/ y/n and that all the eigenfunctions have a similar, sinusoidal, shape near k^. On the 
other hand, the x~^" enhancement of the leading contributions, eq. (I2.24p . is weak because 
the values of u are not large, even for the leading eigenfunctions, and are dropping to ~ 
for larger n, see Fig. [T] We will discuss this more fully at the end of Section 8. 

A reasonable representation of a realistic impact factor of the proton, in terms of the 
eigenfunctions, requires a modification of the assumption of a common phase at some low 
value of transverse momentum. In Section [6l we demonstrate that such phases can be 
obtained with good accuracy by confronting the discrete eigenfunctions with the HERA 
data. 
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3 Dependence of the infrared phase on eigenfunction 
number 

In this section, we explain why the ansatz of a common infrared phase, rj at k = ko is not 
compatible with a "reasonable" form for the proton impact factor such as the one suggested 
in eq. (12. 26 p . We begin by repeating, briefly, the argument in [6] for a constant t], which was 
given within the context of the LO BFKL kernel only. In this approximation the eigenvalue 
equation (with running coupling in the semi-classical approximation) is 

as{k^) f dk'^}Co{k,k')Uk') = cofM- (3.1) 

In the infrared limit, where u <^ Qsik"^) this may be approximated by 

y"rffc'2/Co(k,kO/.(k') = (3.2) 

whose solution is of the form 

Uk) = ^sm{2uo\n{k/ko)+v^), (3.3) 

irrespective of the value of u. Here uq ~ 0.64, is the solution to 

Xo(«^o) = 0, 
Xo being the LO characteristic function. 
A couple of comments are in order: 

1. The argument does not lend itself to an analysis of the kernel which includes the NLO 
contribution in the infrared region where agik"^) becomes large. 

2. For the higher eigenfunctions for which u is indeed negligible compared with the run- 
ning coupling at ko, the phase r] in the solution (13. 3 p can depend on u because the 
NNLO corrections are not taken into account. In addition, non-perturbative effects are 
expected to determine the behaviour of the BFKL amplitude for k < ko, but there is 
no a priori reason to suppose that such an infrared behaviour is given by a common 
phase for all the eigenfunctions — and indeed we show that this cannot be the case. 

The salient feature of the eigenfunctions we obtain is the fact that above the first two 
eigenfunctions the value of fccrit, at which the frequency of oscillation vanishes, is many orders 
of magnitude above where we are fitting data, see Fig. [2] and Fig. [31 This means that in 
the fit region, we are essentially fitting the proton impact factor to oscillatory functions, in 
ln{k/k()), whose frequency actually varies very little as we scan through the relevant range 
of k, as can be seen in Fig. [31 For n > 2, the difference between the eigenfunctions n and 
n + 1 is that the eigenfunction n + 1 has one more oscillation than the eigenfunction n. 
However, in the region of HERA data, k < (9(30) GeV, both eigenfunctions oscillate in a 
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Figure 4: The frequencies z/(/i;) for the eigenfunctions 1, 2, 3, 5, 10, 20, 50 and 120. 

similar way. This is also directly shown in Fig. HI where we plot the frequencies z/(A;) for 
the eigenfunctions 1, 2, 3, 5, 10, 20, 50 and 120. We see that frequencies v vary between 
and 0.75 only. Thus we are effectively performing a Fourier analysis, for which a (generally 
infinite) range of frequencies is necessary, with a very limited frequency range. 

Now, since in this Fourier transform of the impact factor, the quantity v and ln(A;/A;o) 
are conjugate variables, a Fourier transform with a limited range, Az/ ~ 0.75 can only 
be achieved if the function has a correspondingly broad range, Aln(A;/A;o), in the variable 
XYiikjkQ) where 

Az/Aln(A;/A;o) > vr. 

This means that the proton impact factor would have to have non-negligible support up to 



a value /Cmax of transverse momentum, where 



0.75 In 



/co 



~ TT. 



(3.4) 



This gives fcmax ~ 20 GeV, which is much larger than one would want for a proton impact 
factor, which is expected to peak at a few hundred MeV and decrease rapidly thereafter. 

The BFKL kernel possesses a complete set of eigenfunctions and therefore one would 
expect to be able to expand any function (e.g. any proton impact factor) in terms of these 
functions. However, the eigenfunctions which display oscillations with frequency greater 
than the limit discussed above would have negative eigenvalue, lo. Such negative eigenvalue 
solutions are expected to be of negligible importance at sufficiently low x (they contribute 
positive powers of x). They are also very sensitive to unknown higher order corrections to the 
BFKL kernel, as the perturbative expansion in a power series in a^ is strictly only valid in 
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the limit a^ < cu. Furthermore there is no reason to assume that for negative u there exists 
a branch point in u which leads to a discrete spectrum. We therefore seek an expansion 
for the proton impact factor in terms of the (positive u) discrete subset of eigenfunctions, 
which can be achieved provided the boundary conditions of the proton impact factor are 
compatible with the physical properties of this subset of eigenfunctions. The assumption 
of a constant infrared phase rj at some low k = ko imposes a further boundary condition 
on the set of eigenfunctions. The consequent upper limit on the oscillation frequencies of 
the eigenfunctions obtained means that this further boundary condition is not compatible 
with a proton impact factor that decreases on the scale of a few hundred MeV. We are 
therefore led inevitably to the conclusion that we must release entirely the constraint on rj 
and allow this phase to depend on eigenfunction number. This would shift the relative phases 
of oscillations of various eigenfunctions drawn in Fig. H] and allow us to present our impact 
factor in the form of their linear combination. We use a fit to HERA data, to pin down 
the dependence of the infrared phase on the eigenfunction number and thereby understand 
the BFKL amplitude for gluon-gluon scattering for small values of transverse momentum, 
where perturbation theory is no longer valid. 



4 Non- Hermit ian kernel 

The introduction of running coupling into the BFKL kernel, expanded perturbatively as 

/C(k,k') = a,/Co(k,k') + a,'/Ci(k,k') +••• (4.1) 

begs the question as to which transverse momentum (k or k') should be used to determine the 
running. In [H] a so-called "triumvirate" algorithm is employed, which has the advantage of 
being symmetric in k and k' and therefore preserved the Hermiticity of the kernel. Another 
possibility to reach the Hermiticity is to use the similarity transformation [3] 

K{k, k') -^ U{k) K{k, k') U-\k') 

where U{k) depends on a;s(A;^). However, because of the additional dependence of as{k'^) in 
both methods the number of essential eigenfunctions grows substantially, leading to signifi- 
cant technical difficulties. 

In this analysis it was found that it is simpler to work with the set of eigenfunctions 
which are not completely orthogonal 

jd''kU{k)f:{k) = M^n 7^ 5mn. (4.2) 

where, to simplify the notation, we write in the following f^^ = fn- We make the assumption 
that these eigenfunctions form a complete set of states (subject to a given boundary condition 
aX k = ko) so that we may write 

<52(k-k') = 5^/n(k)AC/:,(k'). (4.3) 
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Then the (forward) BFKL equation for the scattering amphtude of a gluon with transverse 
momentum k) off a gluon with transverse momentum k' and centre-of-mass energy ^/s (^ 
/c, /c'), namely 
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^(s,k,k') = 52(k-k')(5(ln^)+ /"A/C(k,q)^(s,q,k') 



(4.4) 



is given by 



^(s,k,k') = $^/™(k)AC;^/„(k')(^ 



(4.5) 




Figure 5: The overlaps Nmm given by eq. (14. 2p . for m = 25. 

As a demonstration, in Fig. [5] we show the overlaps Mmn for m = 25. We can see that 
the off-diagonal elements are much smaller than the diagonal term, as expected from the 
fact that the region where k and k' are substantially different generally contributes little 
to the amplitude. This is a consequence of the fact that in this region the semi-classical 
approximation works fairly well and therefore we can neglect the non-commutativity of 
asik"^) and d/dln{k'^). However, there are many such overlaps and in our analysis we probe 
a region in which k is relatively large, being controlled by the Q^ of the structure functions, 
whereas k' is controlled by the proton impact factor and is expected to be relatively small. 
Thus these off-diagonal components of the amplitude turn out to have a significant effect on 
the quality of the fit obtained. 
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5 Properties of the Green function 



We can define the pomeron Green function in accordance witli eq. (14 .Sp as 

G,(k, k') = J2 /™(k)iV-^/„(k')e-"^ , 



(5.1) 



wfiere y denotes tlie relative gluon rapidity. It satisfies the integral equation 



^^.(k, k') = f £k"Gy> (k, k")Gy-y> (k", k'; 



(5.2) 



for an arbitrary rapidity y' in the interval < y' < y. The quasi-locality and the semi- 
classical approximation ensures that in the integral of eq. (15. 2p . at fixed k, k' and y, the 
essential values of k" are restricted from above. This is a non-trivial property of the Green 
function but not of the eigenf unctions themselves, because in the orthogonality condition 
(14. 2 p large values of k, of the order of fccrit at large n or m, contribute significantly. 

To verify that this quasi-locality property is also present in our numerical evaluation 
we show in Fig. |6] the distribution of the momentum k in the Green function, G'j^(k, k'), 
integrated over k' with the proton impact factor as in eq. (16. 2p . 



Integrated Green function 




k (GeV) 

Figure 6: Distribution of the momentum k in the Green function, ^^^(k, k'), integrated over 
k' with the proton impact factor at y = ln{s/k'^) = ln(l/a; = 10^). 



Figure [6] shows that indeed, when we limit k' < 1 GeV (by integrating with the proton 
impact factor) the contributions of large momenta k are strongly suppressed for /c > 3 GeV. 
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Therefore, the diffusion into large momenta, which could seem natural for large n eigenfunc- 
tions, is strongly suppressed. On the other hand, we will show in subsequent sections that 
large n eigenf unctions make significant contributions and are essential for the description of 
data with high quality. 



6 Comparison with HERA data 

In this section we determine the rj-n relation and the free parameters of the proton impact 
factor by fitting the new, combined, HERA data for the structure function -F2 pLj . In analogy 
to eqs. f l2.23p and f l2.24p we obtain the expression for F2 by performing convolutions of the 
non-Hermitian BFKL amplitude, eq. (14. 5p . with the photon and proton impact factors: 

^if) ^ /f^*p(fc')(^)^'"/^(k'). (6.2) 

The expression for the structure function is then 

F,{X,Q') = J2-^n^-^nMi^^- (6.3) 

m,n 

This expression can be compared directly to data after the rj-n relation is assumed. In the 
search for this relation we were guided by the principle of simplicity and some analogy to 
the Balmer series. In the Balmer series the energies of the quantum levels have a simple 
dependence on the principal quantum numbers n. In the QCD version of the Regge theory 
developed here the BFKL equation is considered to be analogous to the Schrodinger equation 
for the wavefunction of the pomeron. The BFKL kernel corresponds to the Hamiltonian and 
the eigenvalues u to the energy eigenvalues. The BFKL equation determines the dependence 
on u = LOn provided the boundary conditions at small k ~ Aqcd are defined by specifying 
the non-perturbative phase 77 which could depend on n. We have tried several two-parameter 
forms for the n-dependence of rj, but we found that a simple one-parameter functional form 



"max 



gives the best fits to data in terms of x^/^df, where Ndj = A^dat — Apar with A'dat the 
number of fitted data points and A'par the number of fitted parameters. Here, rimax denotes 
the maximum number of eigenfunctions, k is a free parameter and tiq denotes the range of 
rj values used in the fit. The form of this equation guarantees that the phase of the first 
eigenfunction is zero. The value of t]q was determined from the requirement that the value 
of the overlap integral Am should be zero when m = nmax. This condition minimizes the 
contributions of the eigenfunctions with n > nmax, because the frequencies u^k) of the higher 
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n eigenfunctions are very similar in the low k region. We have investigated other functional 
forms for the ra-dependence of rj which include non-zero phases for the leading eigenfunction 
and combinations of both negative and positive powers of n. However, we find that the 
simple form of eq. (16 ■4p gives the best fit. In the fit, only three parameters were determined 
from data: the power k, the normalization constant A and the suppression slope h of the 
proton impact factor, eq. f l2.26p . 

The separate HI and ZEUS inclusive cross-section measurements taken in the period 
1994-2000 have recently been combined to improve accuracy [1]. We fit the 128 measured F2 
data points for neutral-current e'^p scattering with cuts Q^ > 4 GeV^ and x < 0.01. Statis- 
tical and systematic uncertainties are added in quadrature. The overall data normalization 
uncertainty of 0.5% is absorbed into the fitted A parameter in the proton impact factor, 
eq. fl2.26p . We take an input running coupling as{Mz) = 0.1176 with heavy flavour thresh- 
olds at rric = 1.40 GeV, rrih = 4.75 GeV and rrit = 175 GeV. Massive charm and bottom 
quark contributions are included in the photon impact factor [TO] with the same values of 
rric and rUb. 

The computation of large numbers of eigenfunctions together with their overlaps poses a 
considerable numerical problem. Within the adopted computational precision we were able 
to compute up to 150 eigenfunctions and up to 30 overlaps of each eigenfunction with their 
closest neighbours and obtain a value of F2 with a relative precision of ~ 1%. This precision 
is necessary to match the precision of the measured F2 values of ~ 2%. 

The results of the Discrete Pomeron Solution (DPS) flt are shown in Table [^a) as a 
function of the maximal number of eigenfunctions used in the flt, rimax- The flt with only 
the flrst (leading) eigenfunction gives a very bad flt with x^/^df ~ 80. The resulting F2 is 
even qualitatively in disagreement with data. The flt with the flrst flve eigenfunctions starts 
to reproduce qualitatively the main trends of the data although the quality of the flt is still 
poor, x^/^df ~ 3, see Fig. [3 The addition of more eigenfunctions improves the flt quality 
substantially until, with 120 eigenfunctions, a x^/^df ~ 1-2 is achieved. Table [T](a) shows 
an unexpected property of the discrete solution of the BFKL equation: the high quality of a 
flt requires a large number of eigenfunctions. In the third column of Table [T]^a) we show the 
contribution to the x^ from the subset of 43 data points with x < 0.001, where the BFKL 
formalism is most applicable and the flt quality is somewhat better than in the whole region 
of a; < 0.01. In Table [T](b) we show the flt quality if only the 43 data points with x < 0.001 
are included in the flt. In this restricted region, convergence is faster with an increasing 
number of eigenfunctions, as expected, but still at least 40 eigenfunctions are needed for a 
good flt, then convergence is reached after 60 eigenfunctions. The parameters determined in 
the flt for X < 0.001 are similar to the parameters obtained for x < 0.01. However, they are 
less constrained when using only 43 data points (for x < 0.001) compared to 128 data points 
(for X < 0.01), therefore we choose a cut of x < 0.01 for the following results, at the expense 
of requiring more eigenfunctions for an acceptable flt quality. See also the discussion below. 

In Table [2] we show the qualities of flts as a function of the number of overlap integrals, 
^overi, used in the flt. All flts were performed with the same number of available eigenfunc- 
tions, Tijnax = 120. The table shows that the quality of flts improves rapidly when the overlap 
integral correction is included. The flts have only a very small sensitivity to the value of t]q, 
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(a) Fits with cuts of Q^ > 4 GeV^ and x < 0.01: 



^max 


xVN^fix < 0.01) 


xVN^.tix < 0.001) 


K, 


A 


b 


1 


9792 /125 = 78.3 


2123 /43 = 49.4 





156 


30.0 


5 


349.8 /125 = 2.80 


88.8 /43 = 2.07 


3.78 


3.1 ■ 10^ 


78.0 


20 


286.5 /125 = 2.29 


83.3 /43 = 1.94 


0.96 


632 


15.8 


40 


193.3 /125 = 1.55 


54.9 /43 = 1.28 


0.84 


2315 


23.2 


60 


163.3 /125 = 1.31 


44.8 /43 = 1.04 


0.78 


3647 


25.6 


80 


156.5 /125 = 1.25 


43.5 /43 = 1.01 


0.73 


3081 


24.4 


100 


149.1 /125 = 1.19 


41.3 /43 = 0.96 


0.69 


2414 


22.8 


120 


143.7/125 = 1.15 


39.2 /43 = 0.91 


0.66 


2041 


21.8 



(b) Fits with cuts of Q^ > 4 QeV^ and x < 0.001: 



^max 


X^/Ndfix < 0.001) 


K 


A 


b 


1 


894 /40 = 22.4 





193 


29.6 


5 


55.4 /40 = 1.39 


3.65 


4.8 


10^ 


81.5 


20 


52.6 /40 = 1.31 


1.09 


1.1 


10^ 


33.7 


40 


35.5 /40 = 0.89 


1.00 


9.4 


10^ 


68.2 


60 


27.6 /40 = 0.69 


0.91 


1.1 


10^ 


68.2 


80 


27.1 /40 = 0.68 


0.86 


9.3 


10^ 


66.9 


100 


27.3 /40 = 0.68 


0.82 


7.0 


10^ 


64.6 


120 


27.8 /40 = 0.70 


0.79 


6.0 


10^ 


63.6 



Table 1: The quahties of fits, x'^/^df, using up to nmax eigenfunctions with a cut of Q^ > 
4 GeV^ and either (a) x < 0.01 (128 data points) or (b) x < 0.001 (43 data points). In 
the upper table (a) we also show the contribution to x^ given by the 43 data points with 
X < 0.001, without refitting the parameters. The improvement in x'^ obtained by refitting 
the parameters can be seen in the lower table (b). The fits use three parameters: the 
dimensionless power k and the two parameters, A and b, which determine the properties of 
the proton impact factor, and are both given in units of GeV~^. In the photon impact factor 
four flavours are taken into account. The value of rjQ was determined to be rjo = —0.9. The fit 
with nmax = 1 eigenfunction was performed with the optimized constant phase (a) rj = —0.71 
or (b) 1] = —0.69. The fits with one, five and 20 eigenfunctions were performed with one, 
five and 20 overlaps, respectively. All other fits were performed with the maximum number 
of available overlaps, i.e. 30. 
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Figure 7: The results of the fit, performed with one {dotted line) and five {dashed line) 
eigenf unctions, compared to a subsample of the low-x HERA data [1]. 



not more than 1 or 2 units in x^, therefore r^o = —0.9 was used in all fits. 

The fits shown in Table [1] and [2] were made, as in our previous paper [5], using four quark 
fiavours in the photon impact factor. Since the contribution of the bottom quark, although 
small, is present in HERA data we included it in our final fit. The results of the fit performed 
with 120 eigenfunctions, 30 overlaps and five fiavours are shown in Table [3]o 

The final fit achieves x^/^df = 154.7/125 ~ 1.2. This is a very good quality in view of the 
fact that the precision of the data is very high, of the order of 2%. The value of x^/^df ~ 1.2 
means that the precision of the theoretical computation is similar to data precision. This is 
remarkable in view of the fact that NLO corrections to the DIS impact factor are missing. The 
KMS impact factor [10] which we are using takes into account the kinematical constraints, 
which are a part of the NLO correction, but not the complete correction. The lack of full 



^We note that the value of the parameter, b, is such that the proton impact factor eq. (|2.26|) is peaked at 
k ~ Aqcd, as expected. 



^overl 


xVN,f 


K 


A 


h 





354.6 /125 = 2.84 


0.41 


7.80 


1.40 


10 


206.9 /125 = 1.66 


0.50 


69.1 


5.83 


20 


150.8 /125 = 1.21 


0.60 


444.4 


13.5 


30 


143.7 /125 = 1.15 


0.66 


2041 


21.8 



Table 2: The qualities of fits using up to noveri overlap integrals, and the corresponding 
parameters of the fits, with rjQ = —0.9 and four flavours in the photon impact factor. The 
parameters A and b are both given in units of GeV~^. 



xVN,f K A b 


154.7 /125 = 1.24 0.65 1660 20.6 



Table 3: The parameters of the final fit performed with 120 eigenfunctions and 30 overlaps, 
with riQ = —0.9 and five flavours in the photon impact factor. The parameters A and b are 
both given in units of GeV^^. 

NLO corrections to the DIS impact factor could also be responsible for a slight worsening 
of the fit quality when the bottom flavour was added. To support these arguments we also 
performed a fit in a higher Q^ region (Q^ > 8 GeV^, x < 0.01) and obtained a sizeable 
improvement in the quality of the five-flavour (four-flavour) fit, x^/N^f = 1.06 (1.00). The 
parameters of these fits are similar: k, = 0.63 (0.64) and b = 17 A (18.2) GeV^^ for the 
five-flavour (four- flavour) fit. 

We also found that the fit is insensitive to the particular form of the proton impact factor 
as long as we allow the support to be concentrated close to ko. The proton impact factor 
does not explicitly alter the x-dependence as this is always encoded in the eigenvalues of the 
BFKL eigenfunctions. However, the NLO corrections can introduce a scale such that the 
x-dependence associated with eigenfunction n is amended from x"*^" to {x/xq)~'^", thereby 
introducing another free parameter, xq [12]. We have examined this possibility and found 
that the best fits were nevertheless obtained with xq set to unity. 

The best DPS fit with five flavours (full line) is compared to a subsample of data in Fig. |8l 
In the same figure we also show the results of a standard NLO DGLAP global fit (dotted 
line) by MSTW |13j, which included the same combined HERA data [I]. The MSTW fit [13] 
gives a somewhat better x^ = 112 for the 128 low-x HERA data points [T], using around 
30 free parameterqj in a global fit including many other data sets, so the quality of the fit 
to the low-x HERA data is slightly better than the DPS fit. Whereas in our formalism we 
see a considerable improvement in the quality of the fit if we restrict either x < 0.001 or 
Q^ > 8 GeV^, in Ref. [13] the quality of the MSTW fit is similar over such restricted regions 
as it is for the whole region, x < 0.01 and Q^ > A GeV^. The fits to the restricted regions 
give similar fit parameters as for the whole region and the dependence on the number of 



^It is difBcult to judge how many of these ~ 30 parameters in the DGLAP global fit are relevant to the 
low-x HERA data, but a rough estimate is somewhere between five and ten. 
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eigenfunctions, nmax, shows a similar, slow convergence. Therefore we consider the fit to the 
whole HERA low- a; region, for Q^ > A GeV^, as a representative fit and attribute a slight 
worsening of the quality of the fit to other effects which may become significant when x 
is not sufficiently small, such as the effect of the neglected valence-quark contributions, or 
when Q^ is not sufficiently large, such as NLO corrections to the photon impact factor (or 
possible effects from non-linear evolution). 
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Figure 8: The results of the fit, performed with 120 eigenfunctions and 30 overlap integrals 
for each eigenfunction {full line), compared to a subsample of the low-x HERA data [1]. The 
dotted Ime shows the result of the NLO DGLAP global fit by MSTW [13]. 



Figure [9] shows the rate of rise A, defined by F2 oc (1/a;)^ at fixed Q^, as determined in 
the DPS fit and in the direct phenomenological fit to the data [T4]. The present fit describes 
data very well. This is a substantial improvement in comparison to our previous work ^5j, 
where only a qualitative agreement with the observed Q^ dependence of A was achieved. 
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Figure 9: The rate of rise A, defined by F2 oc (1/a;)'^ at fixed Q^, as determined in the DPS 
fit and in the direct phenomenological fit to the data [14j. 

7 r]—uj relation 

The analogy with the Schrodinger equation discussed above suggests that perturbative wave- 
functions can be smoothly extended to low log(fc) values, i.e. into the non-perturbative region. 
In this region an as-yet-unknown dynamics determines the values of the phase of wavefunc- 
tions which in turn determine the boundary conditions rj. The boundaries 1] are determined 
at a low, but still perturbative, value of fco = 0.3 GeVO used in eq. 02.211) . The formal- 
ism described here allows a safe extrapolation of this phase to the non-perturbative region 
k ~ Aqcd by noting that the generalized phase (pi^ik) can be expressed as [6J 



4>uj{k) 



irK(k')\ 



-2i^Jk)\n{k) 



Vt^ik) 



I', 



-1 



{v')dv\ 



(7.1) 



where in the integration by parts the relation //(^(fccrit) = was used. In leading order, the 
integral over dv' is independent of u and k. In NLO this integral becomes u dependent but 
is still independent of fc, because it extends to the highest value of v, which does not depend 
on the value of k^ (see Fig. H]). This allows us to relate the phase rj defined aX k = k^ to the 



^For future studies, it may be wiser to take a larger value of fco and use the infrared oscillation frequency to 
extrapolate the eigenfunctions to smaller values of momentum. This extrapolation is a complicated process 
and so in this paper we have chosen this value of fco as a compromise value. 
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Figure 10: Relation between the non-perturbative phase t^a and the eigenvalues u. The full 
line shows the rj^-uj relation of eq. (I7.4p . 



phase ?7a defined at fc = Aqcd = 0.220 GeV by 



7] = — 2z/(A;o) In 



ko 



A 



QCD 



+ Va, 



(7.2) 



and to determine the rj/^-u relation shown in Fig. [TDl We note that the exact value of Aqcd 
is somewhat arbitrary, depending on the order of perturbation taken for the /3-function and 
the number of active fiavours assumed. However, any change in this value is absorbed by 
a small overall shift in r/A. It does not affect our fit at all since in practice we choose a 
convenient value for ko and use rj rather than r/A- 

It is instructive to express the t]-u dependence shown in Fig. [10] in an analytic form. We 
have shown above that owing to the limited range of the frequency u, near k = koyWe need to 
include a wide range of phases in order to be able to reproduce a proton impact factor which 
is highly suppressed for k ^ Aqcd- On the other hand we know that the eigenvalues, u, 
range from ~ 0.25 for the leading eigenfunction to zero, which is approached as we consider 
more and more eigenfunctions. Since this is a small range of omega, it is inevitable that 
the rj-u relation will contain a singularity as cj — )■ 0. Indeed we note that the approximate 
relation found in our analysis, w„ = 0.5/(1 + 0.95 n), when combined with the rj^n relation 
of eq. (16. 4p , gives 

. 0.0265 , ,„ , 

r/A = 0.4 r-T7-, when u ^ 0. (7.3) 



00 



In addition to these asymptotic terms we found by fitting with a second order polynomial 
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that 



, 0.0265 , , 

i]K = 0.4 - O.Uoj - 1.9lu^ - —^^, (7.4) 



which described the relation between the non-perturbative phase r]\ and the eigenvalues u 
very well, see Fig. [TOl 

In leading-logarithmic approximation (LLA) the du' integral of eq. (17. ip can be evaluated 
analytically as 

where /3q = 11 — 2n//3 and a = f^^' xoi^ldu' ^ 0.92. This led to the original assumption [6J 
of a constant phase, up to nvr as in eq. (12.211) . and gave the relation 

"'^ = :5^ n + , + 1/4- ^'-'^ 



The fact that the RHS of eq. (17.51) diverges as u; — > is a reflection of the fact that 
as the eigenfunction number n — )■ oo, kcnt — >■ oo, and the number of oscillations between 
ko and fccrit also becomes infinite. The simple pole behaviour, however, is only valid in the 
approximation in which the characteristic function, x, is taken only to leading order and the 
coupling is run only in leading order. Once higher order terms are taken into account the 
constant behaviour of the non-perturbative phase t] is converted into a smooth behaviour and 
a power-like singularity (weaker than a perturbative pole). The character of this singularity 
could be sensitive to the value of the /3-function at any point between ko and fccrit- The 
singular term in eq. (I7.4p (with a very small coefficient), obtained from our fit to data, is 
a refiection of the truncation of the perturbative expansion and the possible neglect of any 
BSM physics that may affect the /3-function. 

We recall that in the present evaluation the contributions of all known six quark fiavours 
are taken into account but not the contribution from possible BSM particles such as super- 
partners. The third eigenfunction, which has a critical momentum around 50 TeV, should 
already be sensitive to any such BSM physics. This dependence could be substantial because 
the value of the /3-function is changed significantly in the BSM region and more than 75% 
of the phase integration, eq. (12. 9p . extends over the BSM region, for n > 10 eigenfunctions. 
In other words, the phase at the lowest energy point feg is determined by integration over 
the whole energy region, from ko up to /Ccrit, and since the values of kcnt are very large (see 
Fig. [2]), the phase at ko is mainly determined by the energies beyond BSM thresholds. 

The quantitative effects of BSM physics are very difficult to estimate because they change 
the k dependence of the phases and the running of the coupling simultaneously. Therefore 
any detailed investigation of such effects requires a full NLO evaluation of the BSM effects 
both on the running of the coupling and the BFKL characteristic functions of the BFKL 
equation, which we plan to do in a forthcoming paper. 
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Figure 11: Unintegrated gluon density at x = 10 ^ and x = 10" 



8 Unintegrated gluon density 



The unintegrated gluon density, xg{x, k), is defined in analogy to eq. fl2.24p by 



xg[ 



(^'^) = ^'E(-) fni^)KMi^^- 



U) 



Figure [TT] shows the unintegrated gluon density at two values oi x, x = 10~^ and x = 10^^. 
The unintegrated gluon density determined in the discrete BFKL approach is positive for 
all values of k in the HERA data region. In addition to the results based on fits with 120 
eigenfunctions the figure shows also the results of fits with 40 and 80 eigenfunctions. At x = 
10~^, the fits with 80 (40) eigenfunctions start to differ sizeably from the 120 eigenfunction 
fit at P ~ 100 (10) GeV^. The fit with 80 (40) eigenfunctions has x^ = 156 (193) for 
Ndf = 125 (see Table [1]), which means that the differences from the best fit are very small 
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for 80 eigenfunctions and even for 40 eigenfunctions these differences can be liardly seen 
on the scale of Fig. [HI The fact that differences between the fits with the lower number 
of eigenfunctions can easily be seen in the unintegrated gluon density, at a; = 10~^, but 
not in F2 means that the contribution of higher n eigenfunctions are more important in 
the unintegrated gluon density than in F2. With the deceasing x the convergence with 
increasing number of eigenfunctions improves. At x = 10~^, the unintegrated gluon density 
of the fit with 80 eigenfunctions is almost indistinguishable from the 120 eigenfunction fit. 
The differences between the 40 and 120 eigenfunction fit, although smaller than at x = 10~^, 
are still visible, however. Therefore, HERA data analysed in the discrete BFKL schema 
are determining the unintegrated gluon density up to k"^ < (9(100) GeV^. For higher fc^ one 
should use either more eigenfunctions or change to DGLAP evolution for A;^ > (9(100) GeV^. 

The very good description of the HERA F2 data in the discrete pomeron scheme requires 
a large number of eigenfunctions. This has several reasons. First we recall that the enhance- 
ment factor of the leading pole, x~^'^ , is not very large, e.g. at a; = 10~^ it is only around a 
factor of five, because cji ~ 0.25. Then the contributions of large n eigenfunctions are only 
weakly suppressed in the Green function; the only suppression factor is due to normahzation 
which requires integration over the k region proportional to log(A;crit)- This region grows 
almost linearly with n, which leads to a suppression of the amplitude of the oscillations of 
eigenfunction number n as ~ l/i/n. This, in turn, leads to a slow decrease (also ~ l/^/n) 
of the coefficients of these eigenfunctions required for the fit. Since the proton impact factor 
has a narrow range of support in k it has a large bandwidth and therefore requires a large 
number of eigenfunctions in order to describe with sufficient accuracy the high quality HERA 
data. 

The large number of eigenfunctions required to obtain a good fit to HERA data could 
indicate that the basic QCD interactions which determine the BFKL kernel should be mod- 
ified. Indeed the values of fccrit cross the postulated supersymmetry thresholds (~ 10 TeV) 
between n = 2 and n = 3 and they approach the Planck scale at n ?« 10. This strongly sug- 
gests that the characteristic function at these high energies should be substantially changed 
by crossing thresholds for new physics. The effects of such changes are difficult to estimate 
without performing an evaluation in which the BFKL kernel is modified by N=l to N=4 
supersymmetry effects [T5] . 



At first sight, the suggestion that there can be signals for BSM physics from deep-inelastic 
scattering at energies significantly below the threshold for such new physics appears to be 
in conflict with the general principle of decoupling. However, it is perfectly possible to 
gain information about thresholds for new physics from "low-energy" data coupled with an 
assumption about some high-energy behaviour, generated by the running of couplings. An 
analogy is the prediction of new physics thresholds from the LEP measurements of gauge- 
couplings and the postulate of coupling unification at some GUT scale — these being only 
compatible if the running of the couplings are amended by new physics. In our case, it can be 
argued that the application of the BFKL formalism, which implies a critical momentum for 
which the oscillation frequency vanishes, and the fit at "low-energy", can yield information 
about the running of the QCD coupling between these scales, such that the presence of 
new physics thresholds affects the quality of the fit. It is the infrared phases (universal to 
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all processes mediated by a pomeronjfl that would be affected by such BSM physics and, as 
explained above, such a sensitivity is permitted within the context of the decoupling theorem 
as it is generated by the running of the coupling over a large range. The proton impact factor 
would not display such sensitivity. Since this is unknown, it is also possible that the effect of 
the change in the phases of the BFKL eigenfunctions would be partially compensated by a 
change in the parameters that we use to describe the proton impact factor. Nevertheless, it 
is perfectly plausible that the quality of the overall fit could be sensitive to BSM physics and 
therefore provide some signal, since the shape of the proton impact factor is fairly tightly 
constrained by the requirement of positivity and narrow support. Any detailed investigation 
of these effects requires a full NLO evaluation of BSM effects which we plan to do in a 
forthcoming publication. 



9 Comparison with gluon distribution from a DGLAP 
analysis 

Great care must be exercised if one wishes to attempt to compare the gluon distribution 
obtained here with that obtained from a DGLAP analysis [I6] or even a DGLAP analysis 
supplemented with BFKL dynamics at low x [17]. 



1. In the first place, the comparison can only be made under the assumption that one 
is working at sufficiently low x for the structure functions to be dominated by the 
gluon distribution alone. Although this will eventually be the case if x is sufficiently 
small, in the HERA region the quark distribution used in the DGLAP analyses remains 
numerically important. The gluon contribution is suppressed by an overall factor of 
<^s{Q'^) and also (substantially) by the convolution of the gluon coefficient function 
Cg{x) with the gluon density. The separation of the quark and gluon distributions 
depends on the initial values assumed in the DGLAP analysis at some reference value 
of Q^. In the BFKL approach, it is assumed that there is no primordial quark density 
emerging from the proton and that quarks appear (in NLO BFKL) only from pair 
production from a gluon somewhere along the ladder [3]. This is substantially different 
from the fitted densities extracted from the DGLAP analysis from data used to obtain 
the distributions in P^ITT]. 

2. In the BFKL analysis, the counterpart of the gluon coefficient function is the upper 
impact factor $Dis(a;, k, Q)n The equivalence between this formahsm and the DGLAP 



^At present the only precise data at small enough values of x for the BFKL formalism to be applicable 
are from HERA structure functions. Alternative processes such as low-mass forward Drell-Yan production, 
to be studied at the LHC, may also provide data at the necessary low values of x. 

^In leading order this impact factor is ^-independent, but we use the more careful analysis of |10| in which 
the discrepancy between the measured Bjorken-a; and the longitudinal momentum of the unintegrated gluon 
density is accounted for. 
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approach is obtained within the approximation 

y"^(7(x,P)®$Dis(x,A;,g) ^ Gix,Q^)®Cgix) + as{Q')Pggix)®Gix,Q^)\n(^ 

(9.1) 
where G{x, Q^) is the integrated gluon density, Pqg is an off-diagonal Altarelh-Parisi 
sphtting function, and ® indicates a convolution in x. 

This approximation is obtained in the leading-logarithm approximation, i.e. for suffi- 
ciently large Q^ and the equivalence can be demonstrated in the continuum BFKL case 
(see e.g. [IB]) in the double-leading-logarithm limit, in which the integral over all fre- 
quencies, u, is dominated by the saddle point in the complex u plane. In our analysis, 
we may indeed simulate the continuum case for the higher eigenfunctions, whose eigen- 
values are very closely spaced. However, in order to ensure that the above-mentioned 
saddle point is captured in our discrete sum, we would need to take many more discrete 
eigenfunctions, that number increasing with increasing Q^. This is demonstrated by 
the fact that as we go to higher values of Q^ we need a larger number of eigenfunc- 
tions in order to obtain an integrated gluon density which is stable (convergent) and 
everywhere positive. 

3. Finally, there is the matter of the renormalization prescription dependence of the 
anomalous dimensions themselves, i.e. between the MS usually used in a DGLAP 
analysis and the BFKL prescription to which we are forced in a direct comparison of 
the BFKL analysis with data. This was first considered in [T9] and developed in [20], 
where a translation of the two prescriptions was given which is valid for transverse 
momenta above kcrit but it is unclear how to continue this into the oscillatory region. 
In the analysis of [ITJ this matter was left as an open question. 

We therefore content ourselves with the fact that we have obtained a gluon distribution, 
relevant to our pure BFKL approach, valid up to Q^ = 100 GeV^, which fits HERA data 
very well but which does not lend itself to a comparison with the distributions obtained from 
a DGLAP analysis. 



10 Summary and conclusions 

In this paper we have shown that NLLA solutions of the BFKL equation with the running 
coupling describe all properties of HERA F2 data very well, for Q'^ > 4 GeV^ and x < 10~^, 
provided we allow the infrared phase, 77, to vary with the eigenvalues Un- The solutions of this 
equation have oscillatory form, in which the frequencies i'{k'^) are varying due to the running 
of as{k'^)- We solve the equation near the point z/ = which singles out a specific value of 
k = kcrit, where J^{k'^y.it) = 0. We show that the solutions of the BFKL equation obtained 
here can be considered as a quantized version of the solutions of the DGLAP equation. The 
matching of the BFKL solutions in the region k ~ fccrit leads to a unique set of discrete 
eigenfunctions which cannot be obtained from the DGLAP equation alone. 
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The description of data is obtained by convoluting the Green function with the photon 
and the proton impact factors. The photon impact factor is known, while the shape of the 
proton impact factor was assumed to follow a simple exponential form. The comparison 
with data shows that a particular functional form of the proton impact factor is not very 
important as long as it is positive and concentrated at the values oi k < 0{1) GeV. The 
limited support of the proton impact factor requires, however, a convolution with a large 
number of the eigenfunctions, subjected to a specific phase condition which was determined 
from the fit to data. 

The fitting procedure, especially the finding of the proper infrared phase condition, was 
only possible because of recently published combined HI and ZEUS F2 data from HERA. 
The increased precision of this data requires a large number of eigenfunctions, up to 120, 
to obtain a good fit. The resulting fit permits a very good description of the F2 data and 
the Q^ dependence of the exponent. A, of 1/x, for small-x. The resulting gluon density is 
positive, in the range of HERA energies. 

At higher energies the resulting gluon density is sensitive to the number of eigenfunc- 
tions used in the fit. Since the higher eigenfunctions have eigenvalues which become closer 
together, the inclusion of such eigenfunctions effectively simulates a continuum on top of 
the first few discrete (clearly separated) ones. This could indicate that we are approaching 
a continuum limit which could be better described by the DGLAP evolution alone. So the 
BFKL solution could determine the gluon density up to k"^ of the order of (9(100) GeV^, and 
from then on the DGLAP solution could be used. This could provide a method to overcome 
the problem of negative gluon densities at low x and small scales pertinent to the standard 
DGLAP fits. For example, in Ref. [17], in contrast to the standard DGLAP result, the input 
gluon at Qo = 1 GeV^ obtained from a global fit including small- a; resummation was positive 
and slightly increasing as x — > 0. In view of the importance of the gluon density to the LHC 
physics we plan to study this issue in a forthcoming paper. 

The solutions of the BFKL equation together with HERA data determine the relation 
between the eigenvalue u and the phase r/A which consists of a polynomial term and a singular 
term in u. The polynomial term contains information about the non-perturbative gluonic 
dynamics inside the pomeron because we show that the BFKL equation can be considered 
to be analogous to the Schrodinger equation for the wavefunction of the (interacting) two- 
gluon system. The BFKL kernel corresponds to the Hamiltonian with the eigenvalues Un- 
The analogy with the Schrodinger equation suggests that perturbative wavefunctions can be 
smoothly extended to very low virtuality values, fc^, i.e. into the non-perturbative region. In 
this region an as- yet-unknown dynamics determines the values of the phase of wavefunctions 
which in turn determine the boundary conditions rj/y. 

The singular term, on the other hand, is presumably generated by the perturbative effects 
which were not fully taken into account in our evaluation. This term is sensitive to the high 
virtuality behaviour of the gluon-gluon amplitude, much beyond the virtualities which are 
actually tested in the experiment. This remarkable property is due to the fact that in the 
evolution scheme developed here, the BFKL equation is solved near the critical point, /ccrit, 
and the value of /Ccrit grows quickly with the increase of the eigenfunction number, rapidly 
crossing proposed thresholds for new physics and even the Planck scale. Since we found 
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that the proper description of data requires a large number of eigenfunctions we obtain an 
apparent sensitivity to the BSM effects. We recall that in our approach the eigenfunctions 
with large n correspond to Regge poles which are similar to hadrons with a very small size, 
because fccrit is very large. Their masses (i.e. ujn) are small but can depend on the BSM 
physics. The sheer potential existence of BSM particles, although never produced in the 
interactions relevant to the fitted data, modify the running of the coupling and the (NLO) 
characteristic function of the BFKL equation below the critical momenta and, in turn, modify 
the frequency, amplitude and phase of the eigenfunctions at low-energy. We have shown that 
these states have a soft hadronic tail (i.e. a part of the eigenfunctions around k ~ Aqcd) 
by which they interact with proton and photon and give an essential contribution to the 
structure functions. 

However, only a full NLLA evaluation which takes into account all possible BSM states 
can show whether this apparent sensitivity turns out to be real. If it turns out to be sufficient 
to act as a signal for BSM physics this would provide a new method of "telescoping the Planck 
scale" ED. 
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